Integrand size = 23, antiderivative size = 155 \[ \int \frac {\sqrt {\cos (c+d x)}}{(a+a \sec (c+d x))^3} \, dx=\frac {49 E\left (\left .\frac {1}{2} (c+d x)\right |2\right )}{10 a^3 d}-\frac {13 \operatorname {EllipticF}\left (\frac {1}{2} (c+d x),2\right )}{6 a^3 d}-\frac {\sin (c+d x)}{5 d \sqrt {\cos (c+d x)} (a+a \sec (c+d x))^3}-\frac {8 \sin (c+d x)}{15 a d \sqrt {\cos (c+d x)} (a+a \sec (c+d x))^2}-\frac {13 \sin (c+d x)}{6 d \sqrt {\cos (c+d x)} \left (a^3+a^3 \sec (c+d x)\right )} \]
49/10*(cos(1/2*d*x+1/2*c)^2)^(1/2)/cos(1/2*d*x+1/2*c)*EllipticE(sin(1/2*d* x+1/2*c),2^(1/2))/a^3/d-13/6*(cos(1/2*d*x+1/2*c)^2)^(1/2)/cos(1/2*d*x+1/2* c)*EllipticF(sin(1/2*d*x+1/2*c),2^(1/2))/a^3/d-1/5*sin(d*x+c)/d/(a+a*sec(d *x+c))^3/cos(d*x+c)^(1/2)-8/15*sin(d*x+c)/a/d/(a+a*sec(d*x+c))^2/cos(d*x+c )^(1/2)-13/6*sin(d*x+c)/d/(a^3+a^3*sec(d*x+c))/cos(d*x+c)^(1/2)
Result contains higher order function than in optimal. Order 5 vs. order 4 in optimal.
Time = 1.41 (sec) , antiderivative size = 157, normalized size of antiderivative = 1.01 \[ \int \frac {\sqrt {\cos (c+d x)}}{(a+a \sec (c+d x))^3} \, dx=\frac {\sqrt {\cos (c+d x)} \csc (c+d x) \left (-\left ((241+860 \cos (c+d x)-164 \cos (2 (c+d x))-410 \cos (3 (c+d x))+115 \cos (4 (c+d x))+30 \cos (5 (c+d x))) \csc ^4(c+d x)\right )+520 \operatorname {Hypergeometric2F1}\left (\frac {1}{4},\frac {1}{2},\frac {5}{4},\cos ^2(c+d x)\right ) \sqrt {\sin ^2(c+d x)}+2240 \cos (c+d x) \operatorname {Hypergeometric2F1}\left (\frac {3}{4},\frac {7}{2},\frac {7}{4},\cos ^2(c+d x)\right ) \sqrt {\sin ^2(c+d x)}\right )}{240 a^3 d} \]
(Sqrt[Cos[c + d*x]]*Csc[c + d*x]*(-((241 + 860*Cos[c + d*x] - 164*Cos[2*(c + d*x)] - 410*Cos[3*(c + d*x)] + 115*Cos[4*(c + d*x)] + 30*Cos[5*(c + d*x )])*Csc[c + d*x]^4) + 520*Hypergeometric2F1[1/4, 1/2, 5/4, Cos[c + d*x]^2] *Sqrt[Sin[c + d*x]^2] + 2240*Cos[c + d*x]*Hypergeometric2F1[3/4, 7/2, 7/4, Cos[c + d*x]^2]*Sqrt[Sin[c + d*x]^2]))/(240*a^3*d)
Time = 1.26 (sec) , antiderivative size = 231, normalized size of antiderivative = 1.49, number of steps used = 17, number of rules used = 17, \(\frac {\text {number of rules}}{\text {integrand size}}\) = 0.739, Rules used = {3042, 4752, 3042, 4304, 27, 3042, 4508, 3042, 4508, 27, 3042, 4274, 3042, 4258, 3042, 3119, 3120}
Below are the steps used by Rubi to obtain the solution. The rule number used for the transformation is given above next to the arrow. The rules definitions used are listed below.
\(\displaystyle \int \frac {\sqrt {\cos (c+d x)}}{(a \sec (c+d x)+a)^3} \, dx\) |
\(\Big \downarrow \) 3042 |
\(\displaystyle \int \frac {\sqrt {\sin \left (c+d x+\frac {\pi }{2}\right )}}{\left (a \csc \left (c+d x+\frac {\pi }{2}\right )+a\right )^3}dx\) |
\(\Big \downarrow \) 4752 |
\(\displaystyle \sqrt {\cos (c+d x)} \sqrt {\sec (c+d x)} \int \frac {1}{\sqrt {\sec (c+d x)} (\sec (c+d x) a+a)^3}dx\) |
\(\Big \downarrow \) 3042 |
\(\displaystyle \sqrt {\cos (c+d x)} \sqrt {\sec (c+d x)} \int \frac {1}{\sqrt {\csc \left (c+d x+\frac {\pi }{2}\right )} \left (\csc \left (c+d x+\frac {\pi }{2}\right ) a+a\right )^3}dx\) |
\(\Big \downarrow \) 4304 |
\(\displaystyle \sqrt {\cos (c+d x)} \sqrt {\sec (c+d x)} \left (-\frac {\int -\frac {11 a-5 a \sec (c+d x)}{2 \sqrt {\sec (c+d x)} (\sec (c+d x) a+a)^2}dx}{5 a^2}-\frac {\sin (c+d x) \sqrt {\sec (c+d x)}}{5 d (a \sec (c+d x)+a)^3}\right )\) |
\(\Big \downarrow \) 27 |
\(\displaystyle \sqrt {\cos (c+d x)} \sqrt {\sec (c+d x)} \left (\frac {\int \frac {11 a-5 a \sec (c+d x)}{\sqrt {\sec (c+d x)} (\sec (c+d x) a+a)^2}dx}{10 a^2}-\frac {\sin (c+d x) \sqrt {\sec (c+d x)}}{5 d (a \sec (c+d x)+a)^3}\right )\) |
\(\Big \downarrow \) 3042 |
\(\displaystyle \sqrt {\cos (c+d x)} \sqrt {\sec (c+d x)} \left (\frac {\int \frac {11 a-5 a \csc \left (c+d x+\frac {\pi }{2}\right )}{\sqrt {\csc \left (c+d x+\frac {\pi }{2}\right )} \left (\csc \left (c+d x+\frac {\pi }{2}\right ) a+a\right )^2}dx}{10 a^2}-\frac {\sin (c+d x) \sqrt {\sec (c+d x)}}{5 d (a \sec (c+d x)+a)^3}\right )\) |
\(\Big \downarrow \) 4508 |
\(\displaystyle \sqrt {\cos (c+d x)} \sqrt {\sec (c+d x)} \left (\frac {\frac {\int \frac {41 a^2-24 a^2 \sec (c+d x)}{\sqrt {\sec (c+d x)} (\sec (c+d x) a+a)}dx}{3 a^2}-\frac {16 a \sin (c+d x) \sqrt {\sec (c+d x)}}{3 d (a \sec (c+d x)+a)^2}}{10 a^2}-\frac {\sin (c+d x) \sqrt {\sec (c+d x)}}{5 d (a \sec (c+d x)+a)^3}\right )\) |
\(\Big \downarrow \) 3042 |
\(\displaystyle \sqrt {\cos (c+d x)} \sqrt {\sec (c+d x)} \left (\frac {\frac {\int \frac {41 a^2-24 a^2 \csc \left (c+d x+\frac {\pi }{2}\right )}{\sqrt {\csc \left (c+d x+\frac {\pi }{2}\right )} \left (\csc \left (c+d x+\frac {\pi }{2}\right ) a+a\right )}dx}{3 a^2}-\frac {16 a \sin (c+d x) \sqrt {\sec (c+d x)}}{3 d (a \sec (c+d x)+a)^2}}{10 a^2}-\frac {\sin (c+d x) \sqrt {\sec (c+d x)}}{5 d (a \sec (c+d x)+a)^3}\right )\) |
\(\Big \downarrow \) 4508 |
\(\displaystyle \sqrt {\cos (c+d x)} \sqrt {\sec (c+d x)} \left (\frac {\frac {\frac {\int \frac {147 a^3-65 a^3 \sec (c+d x)}{2 \sqrt {\sec (c+d x)}}dx}{a^2}-\frac {65 a^2 \sin (c+d x) \sqrt {\sec (c+d x)}}{d (a \sec (c+d x)+a)}}{3 a^2}-\frac {16 a \sin (c+d x) \sqrt {\sec (c+d x)}}{3 d (a \sec (c+d x)+a)^2}}{10 a^2}-\frac {\sin (c+d x) \sqrt {\sec (c+d x)}}{5 d (a \sec (c+d x)+a)^3}\right )\) |
\(\Big \downarrow \) 27 |
\(\displaystyle \sqrt {\cos (c+d x)} \sqrt {\sec (c+d x)} \left (\frac {\frac {\frac {\int \frac {147 a^3-65 a^3 \sec (c+d x)}{\sqrt {\sec (c+d x)}}dx}{2 a^2}-\frac {65 a^2 \sin (c+d x) \sqrt {\sec (c+d x)}}{d (a \sec (c+d x)+a)}}{3 a^2}-\frac {16 a \sin (c+d x) \sqrt {\sec (c+d x)}}{3 d (a \sec (c+d x)+a)^2}}{10 a^2}-\frac {\sin (c+d x) \sqrt {\sec (c+d x)}}{5 d (a \sec (c+d x)+a)^3}\right )\) |
\(\Big \downarrow \) 3042 |
\(\displaystyle \sqrt {\cos (c+d x)} \sqrt {\sec (c+d x)} \left (\frac {\frac {\frac {\int \frac {147 a^3-65 a^3 \csc \left (c+d x+\frac {\pi }{2}\right )}{\sqrt {\csc \left (c+d x+\frac {\pi }{2}\right )}}dx}{2 a^2}-\frac {65 a^2 \sin (c+d x) \sqrt {\sec (c+d x)}}{d (a \sec (c+d x)+a)}}{3 a^2}-\frac {16 a \sin (c+d x) \sqrt {\sec (c+d x)}}{3 d (a \sec (c+d x)+a)^2}}{10 a^2}-\frac {\sin (c+d x) \sqrt {\sec (c+d x)}}{5 d (a \sec (c+d x)+a)^3}\right )\) |
\(\Big \downarrow \) 4274 |
\(\displaystyle \sqrt {\cos (c+d x)} \sqrt {\sec (c+d x)} \left (\frac {\frac {\frac {147 a^3 \int \frac {1}{\sqrt {\sec (c+d x)}}dx-65 a^3 \int \sqrt {\sec (c+d x)}dx}{2 a^2}-\frac {65 a^2 \sin (c+d x) \sqrt {\sec (c+d x)}}{d (a \sec (c+d x)+a)}}{3 a^2}-\frac {16 a \sin (c+d x) \sqrt {\sec (c+d x)}}{3 d (a \sec (c+d x)+a)^2}}{10 a^2}-\frac {\sin (c+d x) \sqrt {\sec (c+d x)}}{5 d (a \sec (c+d x)+a)^3}\right )\) |
\(\Big \downarrow \) 3042 |
\(\displaystyle \sqrt {\cos (c+d x)} \sqrt {\sec (c+d x)} \left (\frac {\frac {\frac {147 a^3 \int \frac {1}{\sqrt {\csc \left (c+d x+\frac {\pi }{2}\right )}}dx-65 a^3 \int \sqrt {\csc \left (c+d x+\frac {\pi }{2}\right )}dx}{2 a^2}-\frac {65 a^2 \sin (c+d x) \sqrt {\sec (c+d x)}}{d (a \sec (c+d x)+a)}}{3 a^2}-\frac {16 a \sin (c+d x) \sqrt {\sec (c+d x)}}{3 d (a \sec (c+d x)+a)^2}}{10 a^2}-\frac {\sin (c+d x) \sqrt {\sec (c+d x)}}{5 d (a \sec (c+d x)+a)^3}\right )\) |
\(\Big \downarrow \) 4258 |
\(\displaystyle \sqrt {\cos (c+d x)} \sqrt {\sec (c+d x)} \left (\frac {\frac {\frac {147 a^3 \sqrt {\cos (c+d x)} \sqrt {\sec (c+d x)} \int \sqrt {\cos (c+d x)}dx-65 a^3 \sqrt {\cos (c+d x)} \sqrt {\sec (c+d x)} \int \frac {1}{\sqrt {\cos (c+d x)}}dx}{2 a^2}-\frac {65 a^2 \sin (c+d x) \sqrt {\sec (c+d x)}}{d (a \sec (c+d x)+a)}}{3 a^2}-\frac {16 a \sin (c+d x) \sqrt {\sec (c+d x)}}{3 d (a \sec (c+d x)+a)^2}}{10 a^2}-\frac {\sin (c+d x) \sqrt {\sec (c+d x)}}{5 d (a \sec (c+d x)+a)^3}\right )\) |
\(\Big \downarrow \) 3042 |
\(\displaystyle \sqrt {\cos (c+d x)} \sqrt {\sec (c+d x)} \left (\frac {\frac {\frac {147 a^3 \sqrt {\cos (c+d x)} \sqrt {\sec (c+d x)} \int \sqrt {\sin \left (c+d x+\frac {\pi }{2}\right )}dx-65 a^3 \sqrt {\cos (c+d x)} \sqrt {\sec (c+d x)} \int \frac {1}{\sqrt {\sin \left (c+d x+\frac {\pi }{2}\right )}}dx}{2 a^2}-\frac {65 a^2 \sin (c+d x) \sqrt {\sec (c+d x)}}{d (a \sec (c+d x)+a)}}{3 a^2}-\frac {16 a \sin (c+d x) \sqrt {\sec (c+d x)}}{3 d (a \sec (c+d x)+a)^2}}{10 a^2}-\frac {\sin (c+d x) \sqrt {\sec (c+d x)}}{5 d (a \sec (c+d x)+a)^3}\right )\) |
\(\Big \downarrow \) 3119 |
\(\displaystyle \sqrt {\cos (c+d x)} \sqrt {\sec (c+d x)} \left (\frac {\frac {\frac {\frac {294 a^3 \sqrt {\cos (c+d x)} \sqrt {\sec (c+d x)} E\left (\left .\frac {1}{2} (c+d x)\right |2\right )}{d}-65 a^3 \sqrt {\cos (c+d x)} \sqrt {\sec (c+d x)} \int \frac {1}{\sqrt {\sin \left (c+d x+\frac {\pi }{2}\right )}}dx}{2 a^2}-\frac {65 a^2 \sin (c+d x) \sqrt {\sec (c+d x)}}{d (a \sec (c+d x)+a)}}{3 a^2}-\frac {16 a \sin (c+d x) \sqrt {\sec (c+d x)}}{3 d (a \sec (c+d x)+a)^2}}{10 a^2}-\frac {\sin (c+d x) \sqrt {\sec (c+d x)}}{5 d (a \sec (c+d x)+a)^3}\right )\) |
\(\Big \downarrow \) 3120 |
\(\displaystyle \sqrt {\cos (c+d x)} \sqrt {\sec (c+d x)} \left (\frac {\frac {\frac {\frac {294 a^3 \sqrt {\cos (c+d x)} \sqrt {\sec (c+d x)} E\left (\left .\frac {1}{2} (c+d x)\right |2\right )}{d}-\frac {130 a^3 \sqrt {\cos (c+d x)} \sqrt {\sec (c+d x)} \operatorname {EllipticF}\left (\frac {1}{2} (c+d x),2\right )}{d}}{2 a^2}-\frac {65 a^2 \sin (c+d x) \sqrt {\sec (c+d x)}}{d (a \sec (c+d x)+a)}}{3 a^2}-\frac {16 a \sin (c+d x) \sqrt {\sec (c+d x)}}{3 d (a \sec (c+d x)+a)^2}}{10 a^2}-\frac {\sin (c+d x) \sqrt {\sec (c+d x)}}{5 d (a \sec (c+d x)+a)^3}\right )\) |
Sqrt[Cos[c + d*x]]*Sqrt[Sec[c + d*x]]*(-1/5*(Sqrt[Sec[c + d*x]]*Sin[c + d* x])/(d*(a + a*Sec[c + d*x])^3) + ((-16*a*Sqrt[Sec[c + d*x]]*Sin[c + d*x])/ (3*d*(a + a*Sec[c + d*x])^2) + (((294*a^3*Sqrt[Cos[c + d*x]]*EllipticE[(c + d*x)/2, 2]*Sqrt[Sec[c + d*x]])/d - (130*a^3*Sqrt[Cos[c + d*x]]*EllipticF [(c + d*x)/2, 2]*Sqrt[Sec[c + d*x]])/d)/(2*a^2) - (65*a^2*Sqrt[Sec[c + d*x ]]*Sin[c + d*x])/(d*(a + a*Sec[c + d*x])))/(3*a^2))/(10*a^2))
3.4.91.3.1 Defintions of rubi rules used
Int[(a_)*(Fx_), x_Symbol] :> Simp[a Int[Fx, x], x] /; FreeQ[a, x] && !Ma tchQ[Fx, (b_)*(Gx_) /; FreeQ[b, x]]
Int[Sqrt[sin[(c_.) + (d_.)*(x_)]], x_Symbol] :> Simp[(2/d)*EllipticE[(1/2)* (c - Pi/2 + d*x), 2], x] /; FreeQ[{c, d}, x]
Int[1/Sqrt[sin[(c_.) + (d_.)*(x_)]], x_Symbol] :> Simp[(2/d)*EllipticF[(1/2 )*(c - Pi/2 + d*x), 2], x] /; FreeQ[{c, d}, x]
Int[(csc[(c_.) + (d_.)*(x_)]*(b_.))^(n_), x_Symbol] :> Simp[(b*Csc[c + d*x] )^n*Sin[c + d*x]^n Int[1/Sin[c + d*x]^n, x], x] /; FreeQ[{b, c, d}, x] && EqQ[n^2, 1/4]
Int[(csc[(e_.) + (f_.)*(x_)]*(d_.))^(n_.)*(csc[(e_.) + (f_.)*(x_)]*(b_.) + (a_)), x_Symbol] :> Simp[a Int[(d*Csc[e + f*x])^n, x], x] + Simp[b/d In t[(d*Csc[e + f*x])^(n + 1), x], x] /; FreeQ[{a, b, d, e, f, n}, x]
Int[(csc[(e_.) + (f_.)*(x_)]*(d_.))^(n_)*(csc[(e_.) + (f_.)*(x_)]*(b_.) + ( a_))^(m_), x_Symbol] :> Simp[(-Cot[e + f*x])*(a + b*Csc[e + f*x])^m*((d*Csc [e + f*x])^n/(f*(2*m + 1))), x] + Simp[1/(a^2*(2*m + 1)) Int[(a + b*Csc[e + f*x])^(m + 1)*(d*Csc[e + f*x])^n*(a*(2*m + n + 1) - b*(m + n + 1)*Csc[e + f*x]), x], x] /; FreeQ[{a, b, d, e, f, n}, x] && EqQ[a^2 - b^2, 0] && LtQ [m, -1] && (IntegersQ[2*m, 2*n] || IntegerQ[m])
Int[(csc[(e_.) + (f_.)*(x_)]*(d_.))^(n_)*(csc[(e_.) + (f_.)*(x_)]*(b_.) + ( a_))^(m_)*(csc[(e_.) + (f_.)*(x_)]*(B_.) + (A_)), x_Symbol] :> Simp[(-(A*b - a*B))*Cot[e + f*x]*(a + b*Csc[e + f*x])^m*((d*Csc[e + f*x])^n/(b*f*(2*m + 1))), x] - Simp[1/(a^2*(2*m + 1)) Int[(a + b*Csc[e + f*x])^(m + 1)*(d*Cs c[e + f*x])^n*Simp[b*B*n - a*A*(2*m + n + 1) + (A*b - a*B)*(m + n + 1)*Csc[ e + f*x], x], x], x] /; FreeQ[{a, b, d, e, f, A, B, n}, x] && NeQ[A*b - a*B , 0] && EqQ[a^2 - b^2, 0] && LtQ[m, -2^(-1)] && !GtQ[n, 0]
Int[(u_)*((c_.)*sin[(a_.) + (b_.)*(x_)])^(m_.), x_Symbol] :> Simp[(c*Csc[a + b*x])^m*(c*Sin[a + b*x])^m Int[ActivateTrig[u]/(c*Csc[a + b*x])^m, x], x] /; FreeQ[{a, b, c, m}, x] && !IntegerQ[m] && KnownSecantIntegrandQ[u, x ]
Time = 9.36 (sec) , antiderivative size = 270, normalized size of antiderivative = 1.74
method | result | size |
default | \(\frac {\sqrt {\left (2 \cos \left (\frac {d x}{2}+\frac {c}{2}\right )^{2}-1\right ) \sin \left (\frac {d x}{2}+\frac {c}{2}\right )^{2}}\, \left (348 \cos \left (\frac {d x}{2}+\frac {c}{2}\right )^{8}+130 \sqrt {\frac {1}{2}-\frac {\cos \left (d x +c \right )}{2}}\, \sqrt {-2 \cos \left (\frac {d x}{2}+\frac {c}{2}\right )^{2}+1}\, \operatorname {EllipticF}\left (\cos \left (\frac {d x}{2}+\frac {c}{2}\right ), \sqrt {2}\right ) \cos \left (\frac {d x}{2}+\frac {c}{2}\right )^{5}+294 \sqrt {\frac {1}{2}-\frac {\cos \left (d x +c \right )}{2}}\, \sqrt {-2 \cos \left (\frac {d x}{2}+\frac {c}{2}\right )^{2}+1}\, \cos \left (\frac {d x}{2}+\frac {c}{2}\right )^{5} \operatorname {EllipticE}\left (\cos \left (\frac {d x}{2}+\frac {c}{2}\right ), \sqrt {2}\right )-578 \cos \left (\frac {d x}{2}+\frac {c}{2}\right )^{6}+264 \cos \left (\frac {d x}{2}+\frac {c}{2}\right )^{4}-37 \cos \left (\frac {d x}{2}+\frac {c}{2}\right )^{2}+3\right )}{60 a^{3} \sqrt {-2 \sin \left (\frac {d x}{2}+\frac {c}{2}\right )^{4}+\sin \left (\frac {d x}{2}+\frac {c}{2}\right )^{2}}\, \cos \left (\frac {d x}{2}+\frac {c}{2}\right )^{5} \sin \left (\frac {d x}{2}+\frac {c}{2}\right ) \sqrt {2 \cos \left (\frac {d x}{2}+\frac {c}{2}\right )^{2}-1}\, d}\) | \(270\) |
1/60*((2*cos(1/2*d*x+1/2*c)^2-1)*sin(1/2*d*x+1/2*c)^2)^(1/2)*(348*cos(1/2* d*x+1/2*c)^8+130*(sin(1/2*d*x+1/2*c)^2)^(1/2)*(-2*cos(1/2*d*x+1/2*c)^2+1)^ (1/2)*EllipticF(cos(1/2*d*x+1/2*c),2^(1/2))*cos(1/2*d*x+1/2*c)^5+294*(sin( 1/2*d*x+1/2*c)^2)^(1/2)*(-2*cos(1/2*d*x+1/2*c)^2+1)^(1/2)*cos(1/2*d*x+1/2* c)^5*EllipticE(cos(1/2*d*x+1/2*c),2^(1/2))-578*cos(1/2*d*x+1/2*c)^6+264*co s(1/2*d*x+1/2*c)^4-37*cos(1/2*d*x+1/2*c)^2+3)/a^3/(-2*sin(1/2*d*x+1/2*c)^4 +sin(1/2*d*x+1/2*c)^2)^(1/2)/cos(1/2*d*x+1/2*c)^5/sin(1/2*d*x+1/2*c)/(2*co s(1/2*d*x+1/2*c)^2-1)^(1/2)/d
Result contains higher order function than in optimal. Order 9 vs. order 4.
Time = 0.10 (sec) , antiderivative size = 344, normalized size of antiderivative = 2.22 \[ \int \frac {\sqrt {\cos (c+d x)}}{(a+a \sec (c+d x))^3} \, dx=-\frac {2 \, {\left (87 \, \cos \left (d x + c\right )^{2} + 146 \, \cos \left (d x + c\right ) + 65\right )} \sqrt {\cos \left (d x + c\right )} \sin \left (d x + c\right ) + 65 \, {\left (-i \, \sqrt {2} \cos \left (d x + c\right )^{3} - 3 i \, \sqrt {2} \cos \left (d x + c\right )^{2} - 3 i \, \sqrt {2} \cos \left (d x + c\right ) - i \, \sqrt {2}\right )} {\rm weierstrassPInverse}\left (-4, 0, \cos \left (d x + c\right ) + i \, \sin \left (d x + c\right )\right ) + 65 \, {\left (i \, \sqrt {2} \cos \left (d x + c\right )^{3} + 3 i \, \sqrt {2} \cos \left (d x + c\right )^{2} + 3 i \, \sqrt {2} \cos \left (d x + c\right ) + i \, \sqrt {2}\right )} {\rm weierstrassPInverse}\left (-4, 0, \cos \left (d x + c\right ) - i \, \sin \left (d x + c\right )\right ) + 147 \, {\left (-i \, \sqrt {2} \cos \left (d x + c\right )^{3} - 3 i \, \sqrt {2} \cos \left (d x + c\right )^{2} - 3 i \, \sqrt {2} \cos \left (d x + c\right ) - i \, \sqrt {2}\right )} {\rm weierstrassZeta}\left (-4, 0, {\rm weierstrassPInverse}\left (-4, 0, \cos \left (d x + c\right ) + i \, \sin \left (d x + c\right )\right )\right ) + 147 \, {\left (i \, \sqrt {2} \cos \left (d x + c\right )^{3} + 3 i \, \sqrt {2} \cos \left (d x + c\right )^{2} + 3 i \, \sqrt {2} \cos \left (d x + c\right ) + i \, \sqrt {2}\right )} {\rm weierstrassZeta}\left (-4, 0, {\rm weierstrassPInverse}\left (-4, 0, \cos \left (d x + c\right ) - i \, \sin \left (d x + c\right )\right )\right )}{60 \, {\left (a^{3} d \cos \left (d x + c\right )^{3} + 3 \, a^{3} d \cos \left (d x + c\right )^{2} + 3 \, a^{3} d \cos \left (d x + c\right ) + a^{3} d\right )}} \]
-1/60*(2*(87*cos(d*x + c)^2 + 146*cos(d*x + c) + 65)*sqrt(cos(d*x + c))*si n(d*x + c) + 65*(-I*sqrt(2)*cos(d*x + c)^3 - 3*I*sqrt(2)*cos(d*x + c)^2 - 3*I*sqrt(2)*cos(d*x + c) - I*sqrt(2))*weierstrassPInverse(-4, 0, cos(d*x + c) + I*sin(d*x + c)) + 65*(I*sqrt(2)*cos(d*x + c)^3 + 3*I*sqrt(2)*cos(d*x + c)^2 + 3*I*sqrt(2)*cos(d*x + c) + I*sqrt(2))*weierstrassPInverse(-4, 0, cos(d*x + c) - I*sin(d*x + c)) + 147*(-I*sqrt(2)*cos(d*x + c)^3 - 3*I*sqr t(2)*cos(d*x + c)^2 - 3*I*sqrt(2)*cos(d*x + c) - I*sqrt(2))*weierstrassZet a(-4, 0, weierstrassPInverse(-4, 0, cos(d*x + c) + I*sin(d*x + c))) + 147* (I*sqrt(2)*cos(d*x + c)^3 + 3*I*sqrt(2)*cos(d*x + c)^2 + 3*I*sqrt(2)*cos(d *x + c) + I*sqrt(2))*weierstrassZeta(-4, 0, weierstrassPInverse(-4, 0, cos (d*x + c) - I*sin(d*x + c))))/(a^3*d*cos(d*x + c)^3 + 3*a^3*d*cos(d*x + c) ^2 + 3*a^3*d*cos(d*x + c) + a^3*d)
\[ \int \frac {\sqrt {\cos (c+d x)}}{(a+a \sec (c+d x))^3} \, dx=\frac {\int \frac {\sqrt {\cos {\left (c + d x \right )}}}{\sec ^{3}{\left (c + d x \right )} + 3 \sec ^{2}{\left (c + d x \right )} + 3 \sec {\left (c + d x \right )} + 1}\, dx}{a^{3}} \]
Integral(sqrt(cos(c + d*x))/(sec(c + d*x)**3 + 3*sec(c + d*x)**2 + 3*sec(c + d*x) + 1), x)/a**3
\[ \int \frac {\sqrt {\cos (c+d x)}}{(a+a \sec (c+d x))^3} \, dx=\int { \frac {\sqrt {\cos \left (d x + c\right )}}{{\left (a \sec \left (d x + c\right ) + a\right )}^{3}} \,d x } \]
\[ \int \frac {\sqrt {\cos (c+d x)}}{(a+a \sec (c+d x))^3} \, dx=\int { \frac {\sqrt {\cos \left (d x + c\right )}}{{\left (a \sec \left (d x + c\right ) + a\right )}^{3}} \,d x } \]
Timed out. \[ \int \frac {\sqrt {\cos (c+d x)}}{(a+a \sec (c+d x))^3} \, dx=\int \frac {\sqrt {\cos \left (c+d\,x\right )}}{{\left (a+\frac {a}{\cos \left (c+d\,x\right )}\right )}^3} \,d x \]